1 



Oh! 
(D 

o 



> 
oo 

o 

a^ 
o 
a^ 



X 



Search for Disoriented Chiral Condensates in 158-A GeV Pb+Pb 
Collisions 

Tapan K. Nayak"" * for the WA98 Collaboration 

M.M. Aggarwal"", A. Agnihotri^, Z. Ahammed'^, A.L.S. Angelis'^, V. Antonenko*^, V. Arefiev-^, 
V. Astakhov-^, V. Avdeitchikov^ , T.C. Awes^', P.V.K.S. Baba'^, S.K. Badyal'^, A. Baldine^, 
L. Barabach^, C. Barlag*, S. Bathe*, B. Batiounia-^, T. Bernier^, K.B. Bhalla^ V.S. Bhatia", 
C. BlumeS R. Bock^, E.-M. Bohne*, D. Bucher*, A. Buijs^ E.-J. Buis', H. Biisching*, L. Carlen'" 
V. Chalyshev-^ , S. Chattopadhyay'^, K.E. Chenawi™', R. Cherbatchev*^, T. Chujo", A. Claussen* 
A.C. Das'", M.R Decowski', V. Djordjadze^, P. Donni'', 1. Doubovik^ M.R. Dutta Majumdar'^ 
S. Eliseev°, K. Enosawa", H. Feldmann*, P. Foka'^, S. Fokin% V. Frolov-^, M.S. Ganti^ S. Garpman'", 

0. Gavrishchuk^, F.J.M. Geurts', T.K. Ghoshs, R. Glasow*, S.K. Gupta^ B. Guskov^, H.A. Gustafsson'^ 

H. H. Gutbrod-', R. Higuchi'^, I. Hrivnacova", M. Ippolitov*^, H. Kalechofsky"', R. Kamermans' , 
K.-H. Kampert*, K. Karadjev^ K. Karpio'', S. Kato", S. Kees*, H. Kim^, B.W. Kolb'=, I. Kosarev^, 

1. Koutcheryaev*^, A. Kugler°, P. Kulinich^, V. Kumar^, M. Kurata", K. Kurita", N. Kuzmin-^, 

I. Langbein'^, A. Lebedev"^, Y.Y. Lee'^, H. Lohner 'P , D.P. Mahapatra"*, V. Manko*^, M. Martin'^, 

A. Maximov^, R. Mehdiyev-/', G. Mgebrichvili^ , Y. Miake", D. Mikhalev-'' , G.C. Mishra^ Y. 
Miyamoto", B. Mohanty^, D. Morrison*, D.S. Mukhopadhyay'^, V. Myalkovski-^^ , H. Naef^, 

B. K. Nandi^ S.K. Nayak^ T.K. Nayak^ S. Neumaier*^, A. Nianine^ V. Nikitine-^, S. Nikolaev^ 
S. Nishimura", P. Nomokov-^, J. Nystrand"*, F.E. Obenshain*, A. Oskarsson'", I. Otterlund"* 
M. Pachr", A. Parfenov^, S. Pavliouk^, T. Peitzmann*, V. Petracek", F. PlasiF, M.L. Purschke^ 
B. Raeven', J. Rak°, S. Raniwala'', V.S. Ramamurthy^ , N.K. Rao'*, F. Retiere-^, K. Reygers* 

G. Roland^ L. Rosselef^, I. Roufanov^, J.M. Rubio'', S.S. Sambyal'*, R. Santo*, S. Sato" 

H. Schlagheck*, H.-R. Schmidt'^, G. Shabratova-^ , I. Sibiriak^, T. Siemiarczuk'', B.C. Sinha'^ 
N. Slavine-^, K. Soderstrom™, N. Solomey"', S.P. S0rensen*, P. Stankus^, G. Stefanek'^, P. Steinberg*^, 
E. Stenlund™, D. Stiiken*, M. Sumbera", T. Svensson"*, M.D. Trivedi^ A. Tsvetkov^ C. Twenhofel', 

L. Tykarski"?, J. Urbahn'^, N.v. Eijndhoven', W.H.v. Heeringen', G.J. v. Nieuwenhuizen^ , A. Vinogradov*^ 
Y.P. Viyogi^ A. Vodopianov^, S. Voros'^, M.A. Vos^ B. Wyslouch^ K. Yagi", Y. Yokota", and 
G.R. Young^' 

"'Univ. of Punjab (India) ^Univ. of Rajasthan (India) '^VECC, Calcutta (India) '^Univ. of 
Geneva (Switzerland) '^Kurchatov (Russia) ^ JINR, Dubna (Russia) ^ORNL, Oak Ridge (USA) 
^Univ. of Jammu (India) "^Univ. of Miinster (Germany) ^ SUBATECH, Nantes (France) ^GSI, 
Darmstadt (Germany) ^ NIKHEF, Utrecht (The Netherlands) ^Univ. of Lund (Sweden) '""Univ. 
of Tsukuba (Japan) °NPI, Rez (Czech Rep.) ^KVI, Groningen (The Netherlands) '^INS, Warsaw 
(Poland) ^MIT, Cambridge (USA) ^lOP, Bhubaneswar (India) ^Univ. of Tennessee (USA) 



The formation of Disoriented Chiral Condensates (DCC) [0] in high energy heavy-ion 
colhsions is associated with large event-by-event fluctuation in the ratio of neutral to 
charged pions. The probability distribution of neutral pion fraction, /, in a DCC domain 
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has been shown 0| to follow the relation: 

P(f) = -L_ where f = Ko/N^ (1) 

which is quite different from that of the normal pion production mechanism. The detection 
and study of DCC is expected to provide valuable information about the chiral phase 
transition and vacuum structure of strong interactions. 

Theoretical predictions suggest that the isospin fluctuations caused by DCC would pro- 
duce clusters of coherent pions in phase space forming domains localized in phase space. 
However, experimental studies carried out so far have been restricted to global fluctu- 
ation in the number of charged particles and photons. Here we present first experimental 
results to search for localized domains of DCC based on event-by-event fluctuation in the 
relative number of charged particles and photons detected within the acceptance of the 
WA98 experiment. 

In the WA98 experiment, the main emphasis has been on high precision, simultaneous 
detection of both hadrons and photons. Charged particle hits (Nch) were counted using 
a circular Silicon Pad Multiplicity Detector (SPMD) located 32.8 cm downstream 
of the target with a coverage of 2.35 < r] < 3.75. Photon multiplicities {N^^nke) were 
measured by the preshower Photon Multiplicity Detector (PMD) placed at 21.5 meters 
downstream of the target covering and t] range of 2.9 — 4.2. It consisted of an array 
of plastic scintillator pads, arranged inside 28 box modules, placed behind 3Xo thick 
lead converter plates. A total of 66K central events corresponding to top 5% of the 
minimum bias cross section, determined from the total transverse energy measured by 
the midrapidity calorimeter (MIRAC), were analyzed. Common coverage of PMD and 
SPMD was considered. 

Physics analysis with the cleaned data was performed by comparing with mixed and 
simulated events. Mixed events were constructed from data combining hits from PMD 
and SPMD chosen randomly from different events taking care of the two track resolutions 
while maintaining the correlation between N^^uke and A^^ch, event-wise, as in the data. 

Simulated events were generated using VENUS 4.12 event generator with default pa- 
rameters. The output was then processed through a detector simulation package in the 
GEANT 3.21 framework. This incorporates the full WA98 experimental setup. The effect 
of Landau fluctuations in the energy loss of charged particles in silicon [0] was included 
in the SPMD simulation. For PMD simulation the GEANT results in terms of energy 
deposition in pads were converted to the pad ADC values taking into account the detector 
and readout effects. At this stage the simulated data (henceforth referred to as VENUS) 
and experimental data were ready to be processed with common analysis tools to study 
specific physics issues. 

and N^i, CORRELATION 
The correlation between N^^uke and A^^ch has been studied in smaller 0-segments by 
dividing the 0-space into 2, 4, 8 and 16 bins using the method described in Ref. 0. The 
correlation plots of N^-uke and A'ch have been constructed for each (p bin. A common 
correlation axis (Z) has been obtained by fitting the above distributions with a second 
order polynomial. The closest distance (-Dz) of the data points to the correlation axis 
has been calculated numerically with the convention that Dz is positive for points below 
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Figure 1. (a) Sz distributions for 2,4,8 and 16 divisions in angle. Solid circles are for 
data, solid histograms are for mixed events and dashed histograms are for simulation, (b) 
Corresponding root mean square (rms) deviations. The bars represent systematic errors. 



the Z-axis. The distribution of Dz represents the relative fluctuations of N^-uke and A'^ch 
from the correlation axis at any given bin. In order to compare these fluctuations at 
different scales in the same level, we work with a scaled variable, Sz = Dz/(j{Dz)-, where 
a{Dz) corresponds to VENUS events. 

The Sz distributions are shown in the left panel of Figure 1. The corresponding rms 
deviations for different bins are shown in the right of panel of figure 1. The presence of 
events with DCC domains of a particular size would show up as deviations from general 
correlation line at the proper bin. This would result in a broader distribution of Sz 
compared to those for generic events at that particular division. From figure 1(b) we also 
observe that the rms deviations of data differ significantly from those of mixed events 
when the number of bins in is 8 or less. This indicates the presence of non-statistical 
fluctuations in N^_iike and A^ch for bin sizes > 45*^. 

DISCRETE WAVELET TECHNIQUE 

The novel and powerful method of discrete wavelet transformation (DWT) has also 
been used for the analysis of the data. Simulation studies by Huang et. al. show that 
the DWT analysis could be a powerful technique for the search of localized DCC. For the 
present DWT analysis the azimuthal space of the PMD and the SPMD were divided into 
small bins in 0, the number of bins in a given scale j being The input to the DWT 
analysis is a spectrum of the sample function at the smallest bin in corresponding to 
the highest resolution scale, jmax = 5. Our sample function was chosen to be of the form 

The output of the DWT consists of a set of wavelet or father function coefficients (FFCs) 
at each scale, from j = I,..., {jmax — !)• 

The distribution of FFCs for a generic distribution of particles is gaussian in nature 
^-0. However, the presence of DCC-like fluctuation makes the distribution non-gaussian, 
with a larger rms deviation of the FFC distribution. Comparing the rms deviations of 
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Figure 2. (a) FFC distributions at j = 1 — 4. Solid circles are for data, solid histograms 
are for mixed events and dashed histograms are for simulation, (b) Corresponding root 
mean square (rms) deviations at scales j = 1-4. The bars represent systematic errors. 



the FFC distribution of data, mixed and simulated events one can get an idea about the 
localized fluctuations in the distributions of N^^uke and Nch in the azimuthal space. 

Figure 2(a) shows the FFC distributions at scales j =1-4 for data, mixed events and 
VENUS. The corresponding rms deviations, C,, are shown in figure 2(b). We observe that 
the rms deviations of data, mixed and VENUS events match with each other at scales 
j = 3 and 4. For scales j = 1 and 2, the rms deviations of VENUS and mixed events are 
close to each other whereas there is a clear difference of these values from those of the 
data. In the DWT analysis, the information about the degree of fluctuation at a scale j is 
obtained from FFC distribution at scale j + 1. One can thus deduce that the differences 
seen at j = 2 and 1 are due to the fluctuations in N^-ake and N^t at scale j = 3, i.e. 
around 45° in A0. 

A detailed analysis of the i] — (p phase space distribution of charged particles and 
photons have been carried out using two independent methods. Both of these indicate 
the presence of nonstatistical fluctuations in localized regions in terms of photon and 
charged particle multiplicities. This could have contributions from many different physics 
effects. Attempts are being made to explain this assuming the fluctuations to be solely 
due to the presence of localized domains of DCC. 
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